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Abstract 

Full quantum calculation of Cherenkov gluon radiation by quark and gluon currents and a 
Cherenkov decay of a gluon into a pair of Cherenkov gluons in transparent media is performed. 
Energy losses due to Cherenkov gluon radiation in high energy nuclear collisions are calculated. The 
angular distribution of the energy flow due to the radiation of Cherenkov gluons is analyzed. 

1 Introduction 

Experimental observation of the two-humped structure of dihadron angular correlations in ultrarel- 
ativistic heavy ion collisions at RHIC JJ [2] [3J 2J [5J [3] bearing a remarkable likelihood to the angular 
distribution of Cherenkov photons [7] has brought into the focus of attention a possible existence of the 
phenomenon of Cherenkov radiation of gluons, an idea formulated in [5] [5] and applied to the analysis 
of ring-like structures in cosmic ray events in |10j . 

Interpretation of the experimental data in terms of the Cherenkov radiation of gluons is not unique. 
Theoretical descriptions aiming at describing the double-humped angular pattern of two-particle az- 
imuthal correlations include that in terms of the Mach cone generated by jets propagating in dense 
medium, see e.g. the recent analysis in 11 , as well as in terms of originating from dynamical fluctu- 
ations of the expanding hot and dense fireball \V2\ . From the experimental point of view it has been 
demonstrated [13] that in the case of large rapidity interval between the two particles and for one specific 
choice of transverse momenta bins for trigger and associated particles the resulting angular pattern can 
be completely described by the azimuthal asymmetries of the collective flow. In the case of narrow 
rapidity interval the situation looks different, see a detailed argumentation in [llj . and the problem of 
finding an appropriate description for the experimental data at the level of detalization of [3J is, in our 
opinion, still open. 

The Cherenkov radiation of gluons is a manifestation of nontrivial properties of non-Abelian medium 
created in ultrarelativistic heavy ion collisions 14 . The interpretation of RHIC data in terms of 
Cherenkov gluon radiation and a summary of earlier work was presented in [15) . The analysis of [15] 
was based on a straightforward generalization of the classical Tamm- Frank theory [TB] . A simple field- 
theoretical model of two interacting scalar fields leading to Cherenkov excitations was considered in [T7] . 
A model taking into account the opacity of the medium and rescattering of Cherenkov gluons consid- 
ered in [18] was shown to successfully reproduce the experimental data on double-humped correlations 
[TJ [2J [3J 31 15] . A new line of studies was started in jTHJ [20] where a theory of Cherenkov radiation of 
mesons was constructed in the framework of holographic approach to strong interactions. 

To develop a more reliable theoretical picture for Cherenkov radiation of gluons one has to generalize 
the classical approach of [5J [H [TS] and the simple scalar field model of [T7] to a quantum field theory 
description based on in-medium QCD. The main goal of thepresent paper is to develop such an approach 
to Cherenkov gluon radiation of quark and gluon currentqj. Our consideration is essentially based on 
the quantum theory of electromagnetic Cherenkov radiation developed in [22], see also [23]. Recently 
the approach of [23J was generalized to the case of a moving medium [23] . 

The calculation of Cherenkov gluon radiation by quark currents presented below is a straightforward 
generalization of the Abelian case considered in [223. The calculation of Cherenkov radiation of gluon 
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currents and of the gluon decay into a pair of Cherenkov gluons are new. The corresponding expressions 
and the resulting qualitative picture of the pattern of energy loss related to the Cherenkov radiation 
present the main results of the present paper. 
The plan of the paper is as follows. 

In Section 2.1 we give some general remarks on the physics of Cherenkov radiation. 
In Section 2.2 we compute the rate of the single Cherenkov decay of the quark current. 
In Section 2.3 we compute the rate of the single Cherenkov decay of the gluon current. 
In Section 2.4 we compute the rate of the double Cherenkov decay of the gluon current. 
In Appendix A we describe a simple field-theoretical model justifying the Feynman rules for in- 
medium QCD used in the paper. 



2 Single and double Cherenkov decays 

In this section we compute the spectra of Cherenkov gluons radiated by quark and gluon currents 
and the spectrum of Cherenkov gluons created in the decay of a free gluon into two Cherenkov gluons. 



2.1 General remarks 

The phenomenon of Cherenkov radiation has its origin in the nontrivial changes of the dispersion 
relation for the excitations (quasiparticles) in the medium (as seen from the poles of the propagators): 

(2.1) 



W 2_ k 2 £ ( W ,k) W 2 -k 2 ' 

where s(u>, k) is a (chromo)permittivity of the medium under consideration. In what follows we shall 
concentrate on the simplified treatment in which the standard in-vacuum quark and gluon currents 
interact with the transverse in-medium excitations, the Cherenkov gluonfl In this setting the Cherenkov 
radiation is a decay of a free vacuum particle q(g) into a quasiparticle g and a free particle q(g) possible 
for certain special values of the permittivity e(w,k) > 1 that allow an existence of transverse massless 
excitations, the Cherenkov gluons, so that, e.g., for the Cherenkov radiation of quark current we have 

q(wi, ki) -> q(u 2 , k 2 ) © g{u 3 , k 3 ). (2.2) 

In the simplest QED case the Cherenkov radiation is a decay of a free electron into a free in-medium 
photon and a free electron |22j . Another interesting process to study is a decay of free in-vacuum gluons 

g(ojx, ki) -> g(u 2 , k 2 ) © g(ua, kg). (2-3) 

The Cherenkov gluon emission is of course possible only for special values of energy and momenta 
of the three participating gluons so that the energy-momentum conservation for the considered decay is 
fulfilled. To give a quantitative description for this possibility one has to consider an explicit model for 
the chromopermittivity tensor s(lo, k). Generically chromopermittivity is a nontrivial matrix in the color 
space e ab (w,k). The nontrivial color structure of e a& (w,k) leads, in particular, to the appearance of the 
color Cherenkov rainbow [14]. In what follows we shall confine ourselves to the simplest quasi- Abelian 
case, where £ 0& (w,k) — > S ab e(uj) and use in our qualitative estimates a model for e(w,k): 

e(oj) = e > 1, u) < ljq (2.4) 
e(u) = 1, w > u) - (2.5) 

The Cherenkov radiation is then possible for excitations with energies in the interval uj < luq. 

In what follows we shall use in our numerical estimates the values e = 5 and luq = 3 GeV obtained 
by fitting the experimental data in |18j . The possible physical interpretation of ujq is the border of the 
region of resonances. 



2 Let us note that a more complete treatment of the problem at hand would involve a trilinear interaction of quasiparticles. 
A sketch of the corresponding field-theoretical formalism is given in Appendix A. 
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Figure 1: Cherenkov decay of the quark current q(p) — » q(p — q) + g(q). 



2.2 Cherenkov decay of quark current 



Let us illustrate the approach we use in this paper by presenting a detailed calculation of the spectrum 
of Cherenkov gluons radiated by the massless quark current. The process in question is then a decay of 
a free quark into a free quark and a Cherenkov gluon, 



q{p) -> q(p-q) +g{q), 



(2.6) 



where an incident quark q(p) propagates along the z axis and has the four-momentum p^ = (E, 0, 0, E) 
and g(q) is a Cherenkov gluon having the four-momentum q^ — (uj, |q| sin9, 0, |q| cos (9) and characterized 
by the in-medium dispersion law |q| = yfsuj emitted at the Cherenkov angle 9 with respect to the direction 
of the incident particle. The corresponding cut diagram is shown in Fig. Q[|]. The final quark has the 
four-momentum p /fl = {{E — uj), — (E — uj) sin j3, 0, (E — uj) cos 0). The conservation of four-momentum 
in the decay leads to the following equalities fixing the Cherenkov and recoil angles 9 and /3: 



cos 9 = 



sin/3 



1 



e — 1 uj 
2 E 



e-1 



uj e — 1 
~E ~ 



1/2 



(2.7) 



(2.8) 



E l-u)/E 

The familiar classical expression for the Cherenkov angle cos 9 — 1/y/e follows from (|2.7p in the limit 
lo/E — > 0. Let us note that in the energy range characterizing the trigger and associate particles, 
correspondingly E and uj, in correlation measurements in heavy ion collisions, in particular in reference 
to RHIC data on two-humped azimuthal angular correlations, the energy-dependent term in (|2.7p can 
be numerically important. 

The angle f3 characterizes the straggling of the incident particle in the transverse plane. From (|2.8[) 
we see that in the classical limit uj/ E — > the leading contribution to the transverse momentum of the 
final quark reads 



IP tI 



IP tL/m 



E 



1 

_ \E 



-1 1/2 



(2.9) 



(2.10) 



so that, at given w, the relative importance of transverse struggling is directly dependent on the value 
of e. 

Let us also note that from (|2.7[) there follows the restriction on the energy of Cherenkov gluon 

uj 2 , 
< -= . (2.11) 



E ~ V^+l 

The matrix element for the Cherenkov decay (12.61) q — > qg reads 

U s '(p-q)(-igi l )(t a ) ki u s (p)4 j \q). 



(2.12) 



The polarization vectors of the Cherenkov gluon should satisfy the in-medium transversality condi- 



tion qe^ l \q) = (in the present paper we use the Coulomb gauge) and can be chosen in the forrrQ 



r(l) 



(q) = -^(0,1,0), c 



(2) 



(q) = — p (cos 9, 0, — sin ( 



(2.13) 



3 The actual calculations in the paper are performed by straightforward computation of |M| 2 . 
4 See Appendix A for details. 



Summation and averaging over the spin and color indices of the matrix element squared gives 

2^ E I^W^I 2 = g2( ^ £ ~ 1} (2|p| 2 sin 2 + 2|p||p'|(l-cos/3)), (2.14) 

from which, taking into account the dispersion law for the Cherenkov gluon, it is straightforward to 
compute the differential decay rate into an interval [u),u + duS\: 

As expected, it differs from the QED answer |22j only by the Casimir invariant for the fundamental rep- 
resentation of SU(N C ), Cf = — 1)/2N C , because we chose the simplest model of chromopermittivity, 
which is diagonal in the color space. The corresponding differential energy loss per unit time is simply 
given by 

P^(«|£;)=w7,_>,g(£4E). (2.16) 

The differential energy loss spectrum (|2 . 1 6|) can be used for computing two observables of physical 
interest. 

First, using the fact that there exists a one-to-one correspondence between the Cherenkov angle 9 
and the energy of the Cherenkov gluon u>, it is straightforward to reinterpret (|2.16|) as describing the 
energy flow into an angular interval [8,9 + d0] 

P{6) = P{w{e)[E)^-. (2.17) 

The resulting distribution is shown in Fig. [5J We see that the energy flow is confined to an angular 
interval [9q , 9 C ] (shaded region in Fig. [2]) where the lower limit 9$ is obtained from 

""•-M^'-T^y (2 - is) 

and the upper limit 9 C corresponds to the classical Cherenkov angle cos6> c = l/y/e corresponding to 
taking the limit u/E -> in ([2TT6"]) . 

Second, by integrating the differential spectrum (|2.16[) over lj, one gets an expression for the energy 
loss per unit time: 

^J^ML(E\oj Q ,e) = J dwP q ^ gS (uj\E). (2.19) 

Note that the energy loss per unit time can be easily converted into the energy loss per unit length. The 
loss per unit time and per unit length are connected via the relation 

dE IdE , 

= (2 ' 20) 
where v is the speed of the incident particle. In the chosen system of units the speed of quark and gluon 
is v = 1, so we have the result for the energy loss per unit length 

^21(E[uj ,e) = J dioP q ^ qg {uj\E). (2.21) 

The resulting energy loss is plotted in Fig. [3] We see that the Cherenkov energy loss rate for the 
quark current is quite substantial. 

2.3 Cherenkov decay of gluon current 

Let us now turn to the consideration of the Cherenkov gluon radiation by the gluon current. Analo- 
gously to (|2.6I) this process is a decay of in- vacuum gluon into in-vacuum and Cherenkov gluons: 

g(p) ->ff(p-g)+ g(q), (2.22) 

The corresponding cut diagram for this process is shown in Fig. 2) The kinematics of the gluon 
Cherenkov decay (|2.22[) is completely equivalent to that of (|2.6[) and is described by the equations 

(EMU- 
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Figure 2: The angular differential energy flow 
of quark Cherenkov radiation, e = 5, luq = 
3 GeV and E = 10 GeV. 




Figure 3: The quark Cherenkov energy loss, 
e = 5, loq = 3 GeV. 



The matrix element M. g ^ gg of the decay (|2.22|) reads 



Mtli^ = -igf abc [2(p5 (fc) (q))(e (4) (p)e ( ''' ) (p')) + 



9^99 



+ 2(qeM(p))(e^(p')e W (q))- 



2(qeW(p'))(eW(p)e l " ; (q)) , (2.23) 



(*0, 



where the polarization vectors e ' are described in Eq. (|2.13l) and the conventional in-vacuum polar- 
ization vectors are 



eW(p) = (0,1,0), e( 2 )(p)- (1,0,0), (2.24) 
e (1) (p') = (0,1,0), e (2) (p') = (cos^,0,sin/3). (2.25) 
A straightforward computation leads to the following expression for the differential decay rate 



7s- 



(u>\E) 



a,N, 




-1) 




LJ 


e — 


1 LU 2 \ 




E ~ 


4 


E 2 ) 




+ 1 


H 


E 


)) E2 


£+ l- 


E 


(1 


UJ 

E 





(e + l) 



2 uj 4 



(l-lf E4 



(2.26) 



from which one can compute in complete analogy with the calculations described in the previous para- 
graph. The only difference in the formula for the energy loss 



dE, 



9^99 



(II 



(E\cj ,e) = 



min{u>a,E—uo, } 



dioP, 



9^99 



(2.27) 



where P g ^ gg -{uj\E) 



){lu\E), is the existence of the additional restriction on the energy of the 



emitted Cherenkov gluon uj < E — luq, which is determined by the fact that the energy of the "ordinary" 
gluon after the emission is greater than loq. The angular distribution of the energy flow and the rate of 
the energy loss shown in Figs. [5] and [5] correspondingly. 

The most important feature of the gluonic Cherenkov decay is the large value of the energy loss, see 
Fig. H3 This is to be expected if the cutoff energy u> is not too small and e is not too close to 1 which is 
definitely not the case for the values taken from the fit made in |18j . 

2.4 Double Cherenkov decay of gluon current 

Let us now turn to the analysis of the another purely non-Abelian mechanism coupling ordinary 
gluons to the Cherenkov transverse gluonic excitations, the double Cherenkov decay 



g{p) -> g(p-Q)+9(q), 



(2.28) 
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Figure 4: Cherenkov decay of the gluon current g(p) — > g(p — q) + g{q). 
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Figure 5: The angular differential energy flow 

of gluon Cherenkov radiation, s = 5, uj = Figure 6: The gluon Cherenkov energy loss, s = 

3 GeV and E = 10 GeV. 5, uj = 3 GeV. 



The corresponding cut diagram for this process is shown in Fig. [7] 

Kinematics of the double decay (|2.28[) differs from that of (j2~6j) and (|2.22|) . 

First, the double decay (|2.28[) is possible only fj>r luq < E < 2u>o, i.e. in the restricted interval of the 
energy of the decaying gluon. 

Second, the energy-momentum conservation laws impose restrictions on the decay angles 

g i p 

cos6> = Ve = — (2.29) 

cos/3 = y/i- E , (2.30) 

from which we obtain the following restrictions on the energy of the emitted in-medium gluons 

1 1 lo 1 1 . 
2-27^^2 + VF (2 - 31) 

Let us note that, as follows from (|2.29p . for the typical value E/uj = 2 corresponding to the center 
of the allowed interval (|2.31l) one has cos# = i.e. the angle equals the "classical" Cherenkov one. 

From the restriction that both emitted gluons must be Cherenkov gluons (with energy less than uiq) we 
have the restriction on their energy E — ujq < lo < luq. Then the energy flow is confined to the angular 
interval 

r- £ - 1 E r- E-1E .„„„, 

Jl-—^- ^cos^Kv^-^t^— 2.32 

E - ujq lo 

The boundaries of the angular interval (|2.32[) tend to each other as E approaches 2wo- When E = 2cj , 
the boundaries of the interval coincide, which physically means that the process of double Cherenkov 
decay have no kinematical window to take place. This conclusion is in agreement with the previous 
statement that the double Cherenkov decay takes place only if ujq < E < 2ujq. 
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Figure 7: Double Cherenkov decay of the gluon current g(p) — > g{p — q) + g(q). 
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Figure 8: The angular differential energy flow 
of the gluonic double Cherenkov decay, e — 5, 
uj = 3 GeV and E = 5 GeV. 



Figure 9: The lifetime of the gluon decaying 
through double Cherenkov decay, e = 5, loq = 
3 GeV. 



Calculation of the double Cherenkov decay is completely analogous to that for the Cherenkov decay 
considered in the previous paragraph. The matrix element now reads 

1 

M £$3 = -^f^ [2( P e (fe) (q))(eW(p)e ( ''" ) (p')) + 

+ 2(qe«(p))(5(^(p')e (fc) (q))- 

-2(q^' ) (p'))(e«(p)S (fe) (q))] (2.33) 

and the corresponding expression for the differential decay rate takes the form 



1g^gg{u\E) 



a s N r 



1 E 



E 2 (1 



e l 



1 1 



2e 1 



(2.34) 



The angular distribution of the energy flow for the double Cherenkov decay is shown in Fig. [S] The 
region corresponding to the allowed angular interval ([2. 32ft is shaded. The corresponding energy loss is so 
larg€0 (approximately three times larger than for the single Cherenkov decay considered in the previous 
paragraph) that it makes more sense to plot, instead of it, the lifetime of a decaying gluon as a function 
of its energy shown in Fig. [HI We see that unless the energy is not too close to the decay boundary of 
2ojq, the decay turns out to be extremely fast. 



2.5 Cherenkov energy losses 

The results of studying the single Cherenkov decays of quark and gluon currents and the double 
Cherenkov decay of gluon currents lead to the following generic picture for the Cherenkov-related energy 
loss. 

5 For the initial gluon with the energy E = 5 GeV it is about 35 GeV/fm. 
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• For quark currents the only available decay channel is the single Cherenkov decay. The correspond- 
ing energy loss is non-negligible but subleading with respect to that of the gluon current. 

• For incident gluons with energy in the interval luq < E < 2loq the leading contribution to the energy 
loss comes from the double Cherenkov decay. The corresponding pattern of angular correlations 
corresponds to two peaks around the direction of propagation of the decaying gluon. There also 
exists a small contribution due to single Cherenkov decay. 

• At the threshold energy E = 2uo there takes place a regime switch between the predominant 
double Cherenkov decay at E < 2wo to the single Cherenkov decay of quark and gluon currents 
at E > 2loq where one expects the possible appearance of the third hump corresponding to the 
incident particle. Besides that, as follows from (|2.7p and (|2.29[) . at the threshold E = 2ujo there 
takes place the following change in the Cherenkov angle: 

cos9 = ^-e- £ —l 2 -^. -> cos0 = 4=(l + — — V (2-35) 
2^ u V? V 2 2w o/ 

It is easy to verify that the Cherenkov angle of the double Cherenkov decay is greater than the 
Cherenkov angle for the single Cherenkov decay for all uj kinematically allowed for both processes. 
This means that in addition to the regime switch between the dominant energy loss processes we 
also have a sharp change of the Cherenkov angles of emitted gluons at E = 2wo- 

These features appear to be in qualitative agreement with the pattern of angular correlations mea- 
sured at RHIC [3]. A detailed comparison with the data will be published separately. 
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Appendix 

A Field theory model for in-medium QCD 

In this appendix we consider a simple field theory model justifying the Feynman rules for in-medium 
QCD used in the present paper. Our consideration will be confined to the case of QCD matter in its 
rest system. 

The notion of dielectric permittivity in in-medium QED in the case of the homogeneous, isotropic 
medium with temporal dispersion arises in describing a response to an external electric field: 



/t r+oo 
dt'e(t - t')E(t', x) = / dre(r)E(t - t, x) 
-oo JO 



(A.l) 



or, in the Fourier space: 



r+oo 

D(w,r) = e(w)E(w,r), e(u) = / dre(r)e^ r 

Jo 



(A.2) 



Let introduce the following simple action for the in-medium QCD (note that we use the quasi- Abelian 
model e ab (t - f) = S ab e(t - t')): 



d 4 xTr 



dt'e(t-t') F M (t',x) W(t,t') F 0i {t,x.)W\t,t') - -F^F l] 



(A.3) 



where 



W{t,t') =P\exp 



-igj*A°(T,x)dT 



(A.4) 
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is the Wilson line introduced to preserve the local gauge invariance. It is convenient to choose the 
Coulomb gauge A a0 = 0, d; L A a% — in which the action (|A.3|) simplifies to 



S = I d A x 



\ f dt'e{t - 0<9o^"(i',x)(9o^ ai (£,x) - jF aij F aij 



(A.5) 



Let us now divide the action (IA.5I) into the free field and interaction contributions: 



So 



d 4 x 



d 4 x 



- f dt'e(t - t')d A m {t' ,x)<9 A ai (t,x) - h&A ai - d j A ai ){d l A aj - &>A ai ) 

n 2 

-gf (d i A ai )A A -' — — f abc ^ aed A h% A c 3 A ei A d 3 



(A.6) 
(A.7) 



From the decomposition (|A.6[) there follows that in the Coulomb gauge the form of the triple-gluon 
interaction remains unchanged. Turning now to the free field contribution we obtain after integration 
by parts: 



So = I d^x— 



dt'e{t - t')d A ai (t',x)d A ai {t,x) - (d l A aj )(d i A a °) 



(A. 



It can be easily shown that the expression for the gluon-gluon-Cherenkov gluon and gluon-Cherenkov 
gluon-Cherenkov gluon vertices corresponding to the Feynman graphs is exactly the same as for the 
three gluon vertex in the ordinary non-Abelian gauge theory. Varying the action (|A.8|) with respect to 
the field A al (x) we arrive at the following equations of motion 



or, in the Fourier space, 



dTE(T)(d ) 2 A m (t - r,x) - (a,-) 2 A oi (t,x) = o 



(e(fc°)(fc°) 2 -k 2 ) A ai (k) =0 



(A.9) 



(A.10) 



The above consideration shows that within the chosen model of chromoclcctric permittivity we have two 
different branches of the dispersion relation (two different types of excitations): 



|k| = y/euj, uj < uj 
|k| =u, u > uj 



and, therefore, the following decomposition for A al : 
A a \x) = A$(x) + Afaix) 



A(2) ( x ) 



|k|<^ (2tt) 3 V2F x 

d 3 k 1 
|k|> Wo (2tt) 3 V2F 



-ikx j_ ^af ikx 
^ "k(A) e 



?, A) (k), |k| = -^k": 



k(A) ( 



_L_ ' 

+ C k(A) e 



e\ x) (k), |k| = fc u ; 



(ATI) 
(A.12) 



(A.13) 
(A.14) 

(A.15) 



It is easy to check that in order to be consistent with the commutation relations for the vector potential 
the gluon polarization vectors should be normalized as follows: (e^ 2 ^(k)) 2 = 1 for the "ordinary" gluons 
and s(e(i 2 -,(k)) 2 = 1 for in-medium gluons. 
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